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Topological structure of the SU(3) vacuum and exceptional eigenmodes of 
the improved Wilson-Dirac operator* 
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We present a study of the instanton size and spatial distributions in pure SU(3) gauge theory using under-relaxed 
cooling. We also investigate the low-lying eigenmodes of the (improved) Wilson-Dirac operator, in particular, the 
appearance of zero-modes and their space-time localisation with respect to instantons in the underlying gauge 
field. 
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Instanton content of the SU(3) vacuum ^ 

The importance of the instanton content of 
SU{3) gauge theory comes through both the in- 
trinsic importance of understanding the ground 
state of QCD and the role instantons are conjec- 
tured to play in light hadron structure. Cooling is 
a technique for removing the high-frequency non- 
topological excitations of the gauge-field. How- 
ever, during cooling instantons are also removed; 
either if they are very small (lattice artifacts) or 
through // annihilation. We use under-relaxed 
cooling to reduce the latter problem. Also, on 
the cooled configurations there is still the prob- 
lem of extracting the instanton properties; for 
this we have developed pattern-recognition algo- 
rithms. We present results for 20 configurations 
at /3 = 6.0 (16^48) and (3 = 6.2 (24^48) lattices. 

The gauge update for under-relaxed cooling 
is implemented in each Cabibbo-Marinari sub- 
group as 
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where Umin is the gauge link that minimises the 
action, Uoid the original link, a is the under- 
relaxation parameter and c a normalisation con- 
stant. Under-relaxed cooling increases the num- 
ber of calibrated sweeps needed to annihilate an 
// pair; for a given value of a a calibrated sweep is 



the number of sweeps needed to destroy a p = 2 
instanton. With no under-relaxation one occa- 
sionally finds a very narrow instanton broaden- 
ing out under cooling (presumably because of its 
environment). We have not observed this with 
(significant) under-relaxation. We chose a = 1 
and our measurements were carried out between 
23 and 46 cooling sweeps (corresponding to be- 
tween 10 and 20 cooling sweeps at a = 0). 

On the cooled configurations we first find all 
the local extrema of the symmetrised topologi- 
cal charge density, Q{x), relative to the 3"* block 
surrounding each point. (We do not consider the 
action, S{x), as it clearly records less structure.) 
Each peak is treated as a linear superposition of 
the topological charge of the object at that point, 
calculated from a lattice-corrected formula, plus 
a contribution from every other object on the lat- 
tice, calculated from the continuum formula. A 
self-consistent set of widths is then found by iter- 
ation. These are our candidate instantons. 

Summing up Q{x) over the lattice and compar- 
ing it to nj — Tif shows a discrepancy. We define 



S ^<\Q~ [ni - nj)\ > 
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and impose filters on our candidate instantons. 
The parameters of the filters are chosen to min- 
imise S. We have a "spatial" filter to remove spu- 
rious peaks due to ripples on large objects and a 
"width" filter. The latter compares the width cal- 
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culated above with the width calculated from the 
charge within a radius 2 (or 3) of the peak using a 
lattice-corrected formula; a peak is only included 
if the various widths are in sufficiently good agree- 
ment. Full details will appear elsewhere 

In Figure 1 we show the instanton size distri- 
bution for f3 = 6.2 at 23 sweeps. The distribution 
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Figure 1. Size distribution: /3 — 6.2, 23 sweeps 



is peaked around p w The best fit to the 

large-p tail of the distribution is D{p) cx p^" with 
a ~ 10. As one would expect, the total number 
of instantons is found to vary rapidly with the 
amount of cooling. However the average size and 
the form of the small/large p tails varies much 
less. We note that our results are consistent with 
those of m but not with those of ||]. 

The fact that the impact of a cooling sweep 
does not scale complicates the scaling analysis. 
Figure 2 shows that we can tune the number of 
cooling sweeps so as to get scaling when compar- 
ing P = 6.2 and (3 = 6.0. 

Examining the number of unlike charges a dis- 
tance R away from each peak (normalised by the 
volume of the shell) gives the distribution that 
is shown in Fig 3 (for /? = 6.2 after 23 sweeps). 
It is uniform at long distances and amplified at 
short distances. The corresponding distribution 
for like charges is uniform at long distances and 
suppressed at short distances. This implies some 
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X |S=6.2 46 sweeps - 
o |S=6.0 23 sweeps 
o ^=6.0 28 sweeps - 
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Figure 2. Scaling of size distributions 



Unlike space distribution 
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Figure 3. Unlike Spatial Distribution 



screening of instantons by anti-instantons in the 
vacuum - as expected. 

Calculating < j^^^ > where n(p) is the num- 
ber of objects of size p and q{p) is the charge car- 
ried by objects of size p shows that the charge 
carried by small (large) instantons is correlated 
(anti-correlated) with the sign of Q (Figure 4). 
Indeed our results suggest over-screening of large 
instantons by small anti-instantons. 
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Table 1 

Position kq and chirality x of the zero-mode for 
single instantons without (with) improvement. 
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Figure 4. < ^.gl >, /? = 6.2, 23 sweeps 



Zero modes of the improved Wilson-Dirac 
Operator Q 

The Wilson-Dirac operator, or equivalently 
Q(k) = 75M, may have vanishing or almost-zero 
eigenvalues for certain values of k. Moreover, it 
has been shown that the lowest eigenvalues of Q 
are strongly localised in space-time . 

"Exceptional" configurations are assumed to be 
related to the appearance of (almost-)zero eigen- 
values at some kq < Hcrit- They seem to be 
more frequent at smaller (3 and with SW improve- 
ment. On the other hand, it is unclear if and 
how these zero modes are related to the topol- 
ogy of the underlying gauge field, in particular 
because the chirality x = (V'ItsI^) of the cor- 
responding eigenvectors ^p of Q(ko) is typically 
much smaller than one. To clarify this relation 
we investigate the low- lying eigenmodes of Q(k) 
using a modified conjugate gradient method [Q. 
We verified that the index theorem is realized 
on single instanton configurations (generated as 
in ref. |^), for both improved and unimproved 
Wilson Fermions, with one right-handed (4 right- 
handed plus 3 left-handed) zero modes in the 
range k < 0.2 for anti-periodic (periodic) bound- 
ary conditions. This is in complete analogy to 
the results for 2-dimensional Wilson and for stag- 
gered fermions ^ . 

For anti-periodic boundary conditions, the 
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eigenmodes are localised and centered on the in- 
stanton. Moreover, we find considerable effects 
from discretization errors and SW improvement 
which are summarized in table 1. 

We also investigated the localisation of the low- 
est eigenmodes on four exceptional configurations 
encountered by UKQCD at /? = 6.0 on 16^48 and 
32^64 lattices. In all cases the (almost-)zero mode 
is localized close to a small {p — 2a . . . 3a) instan- 
ton (less than ^/2a away). The chirality has the 
same sign as the topological charge of the topo- 
logical object but is not correlated with the over- 
all topological charge of the configuration. This 
provides some evidence that exceptional configu- 
rations are related to small instantons in the un- 
derlying gauge field. More detailed results will be 
presented elsewhere ||]. 
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